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Abstract 

In this note we report about a method to deal with finite energy sum 
rules. With a reasonable knowledge of the main resonances of the spec¬ 
trum, the method guarantees that we can find a nice duality matching 
between the low energy hadronic data and asymptotic QCD at high en¬ 
ergies. 
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1 Finite Energy Sum Rule 

As a general definition, we can say that a Finite Energy Sum Rule (FESR) is 
an equation that identifies a QCD theoretical calculation along a finite region of 
energy in the physical strong interacting spectrum with the experimental data 
in this corresponding energy region. 

1.1 Two point correlator. 

Our theoretical object is the vacuum expectation value of the time ordered 
product of two currents {j (x)) corresponding to a particular channel (labeled 
by F) of strong interacting particles: 

nr(g) = * J d^xG‘^^{n\T{fix)fio)}\n) (i) 
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This correlator is related to an observable quantity (p^), either scattering 
or decay, through the imaginary part of its analytic structure. The Optical 
Theorem provides the following relationship: 

— Imll^ hadrons) (2) 

TT 

Applying this equation to QCD the first problem is that we cannot perform 
calculations in all ranges of energy. Namely, asymptotic QCD in the low energy 
region where hadron resonances appear, is far from being accessible with the 
present techniques (aside of lattice calculations). Therefore, in order to relate 
low energy properties of resonances (such as masses or decay constants) with 
the high energy QCD parameters (quark masses, strong coupling...) we need to 
introduce appropriate techniques in the correlator to make possible the matching 
between the two regime of energies. A master equation to achieve this target is 
given by the Cauchy’s Theorem of analytic functions. 

1.2 Cauchy’s Theorem 

For an analytic function If'" (s) in the s complex plane with a real cut starting 
in Scut, Cauchy’s Theorem states: 

(f (s) ds — 0 (3) 

27rz / 

C 

C is a closed path surrounding the real cut (figure [^. 



Figure 1: Cauchy’s Theorem on the s plane 
Then, if we introduce an arbitrary entire function K (s) we can also write 
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( 4 ) 


j)K (s) (s) ds = 0 


Being II'" (s) real in the real axes outside the cut we can use Schwarz Re¬ 
flection Principle and write: 
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TT 


rso 

/ K (s) Imll^ (s) ds 

Sriit 


-1 
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K (s) n'" (s) ds 

= Bo 


In this relation we have split the contribution of the integral along the closed 
path C into two parts (figure]^. In the l.h.s. the integration runs along the cut, 
including the low energy region starting at Scuti and in the r.h.s. the integral is 
performed along the circle of radius Sq chosen in the high energy region. Along 
the cut we can use the experimental data information and in the circle we can 
calculate the correlator using asymptotic QCD. In this way we match theory 
with experiment. 

The confrontation of theory and experiment by means of Cauchy’s Theorem 
yields what we call the hnite energy sum rule 


1 

TT 


K (s) Imn^;,,t,,(s)ds 


' ®phys 


1 

27ri 


/|s|=S 0 


^ (■s) IlgcD (■s) 


( 5 ) 


where Sphys stands for the physical threshold. 

Taking into account that the function (®) analytic struc¬ 

ture, with a QCD cut in the real axes starting at sqcd, we can substitute the 
r.h.s. of the former equation by an integral of the imaginary part along the cut, 
thus an equivalent relation is: 


1 1 /'®° 

K{s) ImnJ^^t,^ (s) ds « - / K (s) ImnJ^cD («) ds (6) 

^ d Sphya ^ dSQCO 

Here we have to face the following problem. From the experimental data 
we usually know the physical spectrum near the physical threshold, where reso¬ 
nances appear, whereas at high energies we do not have information of the final 
multiparticle product spectrum. In order to match this information with the 
results of QCD at high energies we introduce a new parameter in the sum rule, 
that we call Ses, above which we are allowed for the substitution of experimental 
data by QCD in the physical cut. 


1.3 The value of s effective (seg) 

Let us define the effective threshold satisfying 

Im nEata (s) ~ ImnJ^cD (s) ^r s ^ Seff 


and 


1 

TT 


pSeii 1 rSeif 

/ a: (s) Im (s) ds «- / 

J Sphys ^ J BQCC 


K isi Tm nlf. 


(7) 

( 8 ) 
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This last relation does not follow from just from the approximation of in¬ 
tegrands Q. We need also to choose suitable “kernels” K (s). As a result we 
face with two problems: how to choose the appropriate kernel and how to hx 
the value for Ses to get a good approximation in equation (|^. 



Figure 2: Physical cut 


2 Choosing the appropriate Kernel 

For sake of comparison we will focus our discussion in two kind of kernels: the 
exponential kernel and the Legendre Polynomial kernels, which are the object 
of this note. 

2.1 Exponential Kernel: SVZ-Laplace Sum Rule 

The exponential kernel, that was first introduced in enhances the resonance 
region in the energy interval [Sphys, Seff] and is given by 

K{s,M^) = 

This falling down of the exponential kernel works in such a way that, on both 
the experimental data and the asymptotic QCD integrals, we can neglect the 
part of the integral in the interval [seffj so] even in the limit sq —>■ oo. 

Then, with this kernel the Sum Rule ([^ becomes 

1 /*Seff 1 fSeif 

- / Im (s) dsK - j Im n{^cD (s, rn) ds (9) 

^ •'Sphya ^ d SQCn 

In the r.h.s. of ^ we have introduced a QCD mass parameter (m) to be 
estimated. The prediction of m is done at some particular M provided that 
there exists a neighbourhood of M where this prediction is constant. Notice 
that M fixes the slope of the exponential kernel in the resonance region, i.e. the 
amount of enhancement or suppression we introduce with the kernel. 
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2.2 Legendre Polynomial (LP) Kernels 

This set of kernels, that was first successfully used in the calculation of the 
heavy quark masses ma, is introduced defining a set of n-degree orthogonal 
polynomials P„(s) in the energy interval s € [seff,so] with a normalization 
condition at the lowest lying resonance, s = M^, as follows: 



Pn{M],) = 1 


q = 0,,l...,n- 1 


( 10 ) 


These polynomials are related to the ordinary Legendre polynomials Vnix) 
in the interval x € [—1,1]: 


Pnis) = 


Vn (xis)) 

Pn{x{Mi)y 


x{s) = 


2s (sq “t“ Seff) 
Sq Seff 


We quote here some of the Legendre Polynomials we use: 


P2{X) 

= ^(3a:'-l), 


PAx) 

= ^(5a:^-3a;), 


Va{x) 

= i(35a:^-30a:^-f 3), 
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Vb{x) 

= i(63a:^ — 70x^ + 15a:). 
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(11) 


Notice that the oscillatory nature of the polynomials cancels the unknown 
data in the interval [ses, sq] more efficiently than in the exponential case. 

We quote here some of the Legendre Polynomials we use: 


^ IsZ. Imn^ata (s) ds 

~ ^ Zqcd ('^5 Seff, So) ImnJ^cD (s. rn) ds (12) 


where (s; Ses, Sq) are defined in the range s S [ses, Sq] according to (10). 

Now the method to determine the QCD parameter m within this sum rule 
works as follows: 


• As a first step we study in equation (12) the stability of our prediction of 
m with So, taking Seg around the continuum physical threshold, beyond 
the resonance region. 


• Then we adjust the value for Seg by demanding optimal stability, namely 
that the function m (so) remains constant in a plateau as large as possible. 

• The estimate for m is taken for the optimal stability values of the param¬ 
eters So and Seg. 
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• Finally, we check the convergence of the result with the degree of the 
polynomial (n). 

To see how this sum rule method works in a particular calculation, we de¬ 
scribe in the next section the determination of the B decay constant (/s). 


3 Legendre Polynomial Kernel in the fs deter¬ 
mination. 


Taking the pseudoscalar correlator for the pair of quarks ub and assuming the 
dominance of the lowest lying resonance {B pole) we take: 

i Im nLta (s) = flM% 5 {s- M%) , (13) 

being the decay constant fs the unknown in equation & 

^ rSeii 

fB (^o) ~ J ®eff) So) (s) ds (14) 


The dependence of /s in Sq is implicit in the definition of the polynomial 
(s; Seff, So) (10). Choosing Seff at the continuum physical threshold, Seff = 
{Mb 


2m,r) = 29.9Gey^ as in IHIl], the polynomials that give stable results 


for different degrees n are depicted in figure 



Figure 3: Legendre polynomials for fs with threshold Sphys = 29.9 Geu^ (see 
equation 0) 

We see that the higher n we choose, the larger values for sq we need to find 
stability, however the slope of the polynomial remains roughly constant in the 
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threshold region. Again we find, as in the exponential case, that the amount 
of enhancement of the resonance region is a fundamental issue in the sum rule 
method. 



Figure 4: fs results for different polynomials at the physical threshold Sphys = 
29.9(see text). 

Also looking at the stability regions for different n degree of the LP kernel we 
appreciate in figure that increasing the degree of the polynomial the plateau 
gets wider and the results, taken at the inflexion points, show a very good 
convergence with the polynomial degree, the result going to the value fs ~ 
175 Mev. 

Next, instead of fixing Seff at the continuum physical threshold, we deter¬ 
mine it by demanding optimal stability [B]. In our case we will achieve this 
by imposing that the first derivative should also vanish at the inflexion point. 
Choosing this way to determine Seg we improve substantially the stability region 
(as it is shown in figure with a result for the decay constant slightly higher 
(/b ~ 186 Mev) which is constant in a wider range of sq- The value of optimal 
stability that we find for Seg is not far from the physical continuum threshold, as 
one could anticipate, although the difference is crucial for stabilizing the result. 
This feature is one of the main motivations to trust the sum rule method with 
LP kernels. 

A further advantage of this method is that the contribution of the sum 
rule integral in the region [seg,so] is really tiny. As a matter of fact we have 
compared the QCD integration pieces in the /b calculation obtaining: 

Pri (s; Seg, 5o) (s, Ul) ds ^ ^ 

/sQCD (S: ds 

The same would happen in the experimental data integration provided we take 
Seg beyond the resonance region. Obviously the lack of precise information 
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Figure 5: Jb results for different polynomials at optimal Se// = 31.3 Gei;^ (see 
text). 


introduces a systematic uncertainty of the method which is beyond our control. 
Nevertheless, we expect this uncertainty to be small, since even the resonances 
in the continuum region are substantially suppressed by the polynomials with 
respect to the lowest lying resonance (see figure . 

4 Legendre Polynomial Kernel in the determi¬ 
nation of the strange quark mass 

The calculation of the light quark masses from the pseudoscalar current is a bit 
cumbersome due to the poor convergence of the QCD correlator with respect to 
the strong coupling. We give here some preliminary results using the LP sum 
rule method. 

For the strange quark mass we use in the spectral function, aside from the 
Kaon pole, the contribution of the K(1460) and K(1830) resonances that sig¬ 
nificantly improves the convergence of the results. In order to perform the 
integration of the experimental side we use a Breit-Wigner model, as depicted 
in figure 

To have a flavour of the results obtained we present in figure the stable 
values for the quark masses obtained with the 6th order Legendre polynomial. 
The stability regions of the mass with the value of sq is apparent, presenting 
again a nice plateau around the stability points. Results for different orders in 
the strong coupling constant show a fairly good convergence of the results for 











Figure 6: Resonance model in the Kaon channel (the kaon pole is omitted). 


the masses. This convergence is worse than in the light-heavy quarks system 
calculations and further investigation is under way. 

Aside from this consideration, what is important to stress here concerning 
the LP sum rule method, is that the value for Ses is determined by optimal 
stability and it is located, as expected, just after the resonance region where 
one expects that the experimental results approach to some smooth function of 
the energy. On the other hand we have checked that the result has a very good 
convergence when the degree of the polynomial is increased. 

5 Conclusions 

To summarize, we have reviewed the main features of the FESR method with 
LP kernels. The main advantages that we find are: 

• The LP kernels eliminate very efficiently the contribution of asymptotic 
QCD and the experimental data in the interval [seff, sq] of the Sum Rule 
integral. 

• The LP kernels are easy to integrate with asymptotic QCD in the circle 
|s| = So and one does not need to extract the imaginary part from the 
QCD correlator. 

• The method is able to determine Sgff in a systematic way, which provides 
a better stability of the results (figure . 

• The final results show nice plateaus with sq as well as a good convergence 
with the degree of the polynomial. 

• The slope of the LP near the threshold approaches to a constant when 
increasing the degree of the polynomial. 
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• First used in the calculation of the heavy quark masses Elia we have used 
it recently in the calculation of the meson decay constants in the heavy- 
light quark systems [aElle] and work in the pion and kaon channels is in 
progress |7]- 
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Figure 7: TOs + rriu results at different orders of the strong coupling constant 
with a 5th degree polynomial {sgff = 3.6 Geu^). 
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